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ABSTRACT 


This  report  traces  the  development  of  a  differential 
correction  theory  for  lunar  and  space  probes.  To  overcome  the 
singularities  Inherent  In  differential  processes  with  conventional 
two-body  descriptions  of  the  pith,  as  the  eccentricity  passes 
through  unity  In  either  the  parabolic  or  rectilinear  sense,  a  uni¬ 
fied  formulation  of  the  two-body  equations  was  developed.  The 
differential  correction  procedure  Is  based  upon  an  ephemerls,  Inte- 
, grated  In  Encke  form  by  special  perturbations,  and  upon  differential 
expressions  derived  from  the  unified  two-body  equations.  The  report 
Includes  experimentation  with' simulated  range,  range-rate,  azimuth, 
altitude,  right  ascension  and/or  declination  data  to  evaluate  the 
performance  of  the  resulting  lAl  709  program. 
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SECTION  1 
INTRODUCTION 


Under  an  amendment  to  contract  AF19 (604) -5885, 
Aeronutronic  has  undertaken  the  development  of  a  differential 
correction  program  for  geocentric  orbits  of  medium  to  high 
eccentricity,  including  lunar  probes.  Orbits  of  interest  in¬ 
clude  the  parabolic,  elliptic,  hyperbolic,  and  rectilinear 
forms,  and  an  orbit  theory  which  was  applicable  at  once  to  all 
of  these  forms,  without  the  singularities  normally  associated 
with  the  transition  between  these,  was  sought.  This  report 
discloses  the  development  of  a  unified  orbit  theory  especially 
adapted  to  machine  computation)  included  are  the  necessary 
differential  relationships  and  details  of  the  Encke  perturbation 
method . * 


*  The  proposal  to  employ  a  unified  formulation  of  the  two-body 
equations,  and  the  major  portion  of  the  development  thereof, 
are  due  to  Samuel  Herrick.  Others  on  the  staff  of  Aeronutronic 
who  participated  in  this  development  include  Jeannine  Arsenault, 

R.  H.  Gersten,  G.  Matlln,  C.  Tross,  C.  T.  Van  Sant,  G.  6.  Westrom, 
and  L.  G.  Walters. 
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SECTION  2 


TECHNICAL  APPROACH 


The  Keplerlan  or  two-body  orbit  is  the  basis  for  the 
perturbation  methods  of  orbit  integration  and  for  the  differential 
expressions  required  for  orbit  correction.  In  the  perturbation 
method;  the  Keplerian  path  is  employed  as  a  reference  path,  and 
only  the  non  two-body  accelerations  (or  "perturbations")  need  be 
integrated  to  define  the  motion.  Differential  expressions  involve 
the  determination  of  linear  cause-and-ef f ect  relationships  between 
observation  and  orbit;  and  require  a  description  of  the  orbit  free 
of  singularities  in  the  region  of  interest.  In  both  the  perturbation 
method  and  differential  correction  use  of  other  than  two-body  (or 
nearly  two-body)  form,  such  as  position  and  velocity,  severely  limits 
the  computational  efficiency. 

A  review  of  the  two-body  orbit  description  betrays  two 
serious  singularities  in  the  region  of  interest.  First  is  the 
transition  from  ellipse  to  parabola  to  hyperbola  as  the  eccentricity 
passes  through  unity  and  the  semi-major  axis  becomes  infinite.  This 
region  is  described  by  three  different  sets  of  formulae,  none  of 
which  is  valid  in  this  entire  region.  The  situation  is  similar  to 
the  engineer' s  description  of  damped  harmonic  motion;  wherein  three 
mathematical  descriptions  cover  the  undamped,  critically  damped, 
and  overdamped  cases.  A  more  subtle  singularity  occurs  when  the 
eccentricity  passes  through  unity  but  the  semi-major  axis  remains 
finite.  The  two-body  orbits  are  degenerate  forms  of  the  ellipse 
and  hyperbola,  but  the  path  is  confined  to  a  rectilinear  path  and 
a  description  which  Involves  the  orientation  of  the  orbit  plane  fails. 
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The  most  familiar  description  of  two-body  motion  is  in 
terms  of  the  six  elements: 

a  semi-major  axis 

e  eccentricity 

T  time  of  perigee  passage 

i  Inclination 

^  longitude  of  the  ascending  node 

^  argument  of  perigee 

This  set  of  elements  embodies  the  following  singularities; 

1.  For  zero  eccentricity  (i.e.,  the  circular  orbit) 
cuis  indeterminate^  and  the  anomalies  used  to  describe  position 

in  the  orbit  plane  have  an  indeterminate  origin. 

2.  For  zero  inclination  (i.e.,  the  equatorial  orbit),  the 
longitude  of  the  node  is  indeterminate,  and  the  argument  of  perigee 
has  an  indeterminate  origin. 

3.  For  rectilinear  motion,  the  orientation  of  the  orbit 
plane  is  undefined)  consequently,  both  the  longitude  of  the  node  and 
the  inclination  are  indeterminate. 

4.  For  nearly  parabolic  motion,  the  equations  for  ellipse 
and  hyperbola  are  singular  for  unit  eccentricity,  whereas  the  para¬ 
bolic  equations  are  invalid  for  non-unit  eccentricity. 

For  the  latter  two  reasons  a  description  of  the  two-body  motion  by 
'these  familiar  elements  must  be  abandoned,  and  an  alternate  set  of 
parameters  chosen  for  this  study.  Incidentally,  Aeronutronic  has 
developed  descriptions  which  also  overcome  singularities  of  the  first 
two  types)  results  of  this  research  are  published  elsewhere.* 


See,  for  example,  Aeronutronic  publication  U-880. 
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The  development  of  the  orbit  theory  is  given  in  the 
following  sections.  These  include  the  basis  for  parameter  selection, 
the  Encke  ephemeris  formulation  in  terms  of  these  parameters,  the 
differential  expressions  for  orbit  correction,  and  the  geometric 
transformations  required  for  the  utilization  of  topocentric  obser¬ 
vation  residuals. 


2.1  SELECTION  OF  PARAMETERS 


Singularities  of  the  types  noted  above  arise  from  the 
parameters  employed  to  describe  the  two-body  orbit.  This  section 
will  introduce  a  set  of  six  parameters  which  removes  all  singularities 
save  that  for  zero  eccentricity  from  the  theory,  thereby  unifying 
the  separate  descriptions  heretofore  required  in  the  vicinity  of  unity 
eccentricity*. 


To  provide  continuity  from  the  elliptic  to  the  hyperbolic 
region.  Including  the  parabolic  (e  s  1,  infinite  a)  and  rectilinear 
(e  >  1,  finite  a)  cases,  it  is  necessary  to  adopt  a  uniform  concept 
of  mean  motion  n,  ordinarily  defined  by; 


"e  '  nTm 

np  =  k'^^M 
nj^  *  k-Ny/i 


elliptic  case 
parabolic  case 
hyperbolic  case 


In  the  following,  the  parabolic  definition  for  mean  dally  motion  will 
be  adopted  for  all  three  cases,  and  the  remaining  two -body  equations 
will  be  adjusted  accordingly.  From  this  procedure  evolves  a  set  of 
two-body  equations  Identical  for  the  elliptic,  parabolic,  rectilinear, 
and  hyperbolic  cases,  with  no  singularities  as  the  eccentricity  passes 
through  unity. 


Parallel  research  undertaken  by  Dr.  Samuel  Herrick  of  Aeronutronic 
has  lead  to  an  alternate  development  in  terms  of  closed  form  f  and  g 
series  expressions  which  is  valid  for  all  eccentricities,  including  zero, 
and  all  inclinations.  Further  development  of  this  completely  unified 
theory  is  continuing;  the  development  of  both  a  variatlon-of-parameters 
theory  and  of  differential  expressions  for  orbit  correction  will  be  pursued. 
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Beglnnlijg  with  the  parabolic  definition  of  mean  dally 
motion^  designated  n^*  the  mean  anomaly  Is  related  to  time  by 


M  -  M  =  n  (t-t  ) 
o  '  o 

-3/2  ~ 

where  the  usual  elliptic  definition  of  M  Is  a  M.  Kepler's  equation^ 

In  elliptic  form^  Is  written 


M  ■  E  -  e  sin  E  «  «  (E  -  sin  E)  +  (1-e)  E 


or 


M  -  M  =  e  U  +  q  X 

where^  by  definition^ 


A 

X  « 

•f~a'E 

u  > 

A 

31 

'^5 

X 


'^7 

X 


31  a  .  5'.  a  •  7! 


The  function  U  is  derived  from  the  (E-sin  E)  series)  for  negative  a 
(hyperbolic  case)  this  series  ma^  be  derived  from  (slnh  F  -  F)  series. 
For  Infinite  a,  (parabola),  the  M  form  for  Kepler's  equation  reduces 
to  the  familiar  Barker's  equation  for  the  parabola.  Thus  the  single 
form  Is  adequate  for  all  conics. 


In  this  unified  development,  the  tilde  (  '^  )  symbol  will  be  used 
to  denote  quantities  whose  definitions  differ  by  virtue  of  the  unification 
process,  from  familiar  quantities  in  the  conventional  two-body  theory. 

This  comment  does  not  apply  to  the  usage  In  Section  2.4,  where  the  tilde 
symbol  Is  used  to  distinguish  between  equatorial  and  horizon  coordinates 
of  the  unit  vector  triad  L,  A  and  D. 
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The  position  vector  £  Is  readily  related  to  the  above 
quantities  through 


I  ‘  ^  I  +  yuj  S- 

To  avoid  the  Indeterminate  nature  of  Q  for  the  rectilinear  case,  where 
the  seml-latus  rectum  p  ■  a  (1-e^)  Is  eero  and  the  orbit  degenerates 
to  a  single  line  along  P,  a  vector  ^  avp  la  introduced  such  that 

-  "  +  -z —  fl- 

V  p 

The  coefficients  are 


r  a  (cos  E  -  e)  X  a  (1-e)  -  a  (1-cos  E) 


=  q  -  C 

where  C  ■  -  -  -  +  — s—  -  . , . 

21  a  •  4 I  a  • 6 ! 


^0)  _  1  ~ 
- —  -VaslnExX  -  -  U 


Note  that  the  elliptic  expressions  above,  expressed  In  terms  of 
circular  functions,  become  hyperbolic  functions  for  negative  a, 
typical  of  the  formulae  for  the  hyperbola.  The  velocity  expression 
r  Is 
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where 


^  - 


=  -  ^  sin  E  = 
r 


r  a 


k.'S.  CO.  E=^  (1  -ic) 

r  r  ^  a  ^ 


The  unified  theory  presented  here  lends  itself  to  automatic  machine 
calculation,  where  transcendental  functions  are  routinely  computed 
from  their  ^aylor  series  representations;  on  the  other  hand,  lack  of 
tables  for  U  and  7!  in  terms  of  and  a"^  will  discourage  hand  computa 
tlon . 

The  parameters  which  suggest  themselves  for  representation 
of  all  orbits  (save  zero  eccentricity;  note  the  usage  of  P)  are 


Mq.  a  .  £.  a  • 

The  previous  development  indicates  the  procedure  for  translating  from 
parameters  to  position  and  velocity.  The  inverse  transformation  from 
position  and  velocity  to  these  parameters  follows: 

•  ^  ^ 

Given  r,  r,  and  t,  to  obtain  M^,  ~,  P  and  2^ 

r^  =  r  •  r 

rr  =  r  •  r 
•2 

s^  =  r  •  r 


Compute  —  and  q; 


1  2 

~  ~  (vis  -  viva  Integral) 

p  =  (r2  /jU 

e^  =  1  -  p/a 
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q  =  p/(l  +  e) 


Compute  P  and 


\ 

I 

Compute  and  ft":* 


~  1 

C  =  J  (r  -  q) 
S  =  rr/e'/F 


S-tiSl  r,+  tjt  (,.»)]  i 


"5  =  if- 


a. 3! 


V 


r\j  rxj 

M  =  q  X  +  e  U 


Compute  M^: 


n  = 

/\/  ^  r\j 

Mq  =  M  -  n  (t  -  tp) 


2.2  ENCKE  EPHEMERIS  PROGRAM 

The  numerical  integration  of  the  perturbative  accelerations, 
(and  of  the  total  acceleration  in  Cowell's  method),  is  designated  as 
"special  perturbations".  Factors  to  be  considered  in  an  orbit  theory 
based  on  special  perturbations  include:  what  is  to  be  integrated  numer¬ 
ically,  the  choice  of  reference  orbit,  and  the  choice  of  the  independent 
variable.  In  addition,  the  type  of  numerical  integration  procedure, 
e.  ,Adams-Bashforth  second-difference,  second-sum  Runge-Kutta,  etc., 
must  also  be  considered. 

The  special  perturbations  approach  is  particularly  well 
suited  to  space  trajectories  which  depart  only  slightly  from  orbits 
for  which  an  analytical  solution  exists,  e.g.,  a  lunar  trajectory. 


rV 

*  Alternatively,  series  for  U  directly  in  terms  of 
be  employed,  i.e.. 


S  may 


and 


^ 

3.'  a-5.’ 


a^‘7.' 


/V  A/  ^ 

M  »  U  +  qS 
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In  connection  with  what  is  to  be  integrated,  a  variety  of 
methods  for  special  perturbations  have  been  developed.  In  effect, 
however,  all  of  these  methods  use  some  modification  or  combination  of 
the  basic  concepts  characterized  by  (1)  the  integration  of  the  totals 
acceleration,  as  in  Cowell's  method,  (2)  integration  of  the  departures 
from  a  fixed  reference  orbit,  as  in  Encke ' s  method,  and  (3)  integration 
of  the  parameters  of  a  varying  reference  orbit. 

The  Encke  ephemeris  program  has  been  shown  to  be  highly 
efficient  for  lunar  trajectories.  The  program  has  demonstrated  its 
ability  to  integrate  trajectories  to  the  moon  with  15  integration 
steps  with  less  than  five  miles  accumulated  truncation  error  at  impact. 
Similar  calculation  with  Cowell's  method  required  almost  200  integra¬ 
tion  steps. 

In  the  Encke  orbit  determination  method  only  the  departures 
from  a  two-body  reference  orbit  are  integrated.  The  position  at  some 
epoch  (t)  on  the  reference  orbit  is  given  by  r  .  The  deviatiorfsof  the 
vehicle  from  the  reference  orbit  can  then  be  given  as  (see  Fig.l): 

jO  -  r  -  r  2.2.1 

~e . 

Considering  one  component  we  have  the  derivatives; 


where 


Finally 


d  s  d  X  d  X 

d'^  dT^  dr^ 


•.  •> 

X  +  X 


-  _Mx 

3 

r 

perturbations 


2.  2.  2 
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FIG.  1 


DEVIATION  FROM  REFERENCE  ORBIT 


1 


For  small  departures  from  the  reference  orbit,  as  in  the  case  of  lunar 
trajectories,  the  first  term  on  the  right  side  of  the  above  equation 
is  small  and  hence  one  integrates  only  slightly  more  than  the  x 
term.  The  first  term  on  the  right  side  is  normally  developed  into  a’ 
series  to  avoid  obtaining  small  differences  of  large  quantities. 

This  series  is  developed  in  the  following  manner; 

First,  X  -  X  +  etc.,  are  substituted  into  the  defining 
equation  for  r2 

2  2  2  2 
r  «  X  +  y  +  z  , 

yielding 

-  (Xg  +  ^  )^  +  (y^  +  Tj)^  +  (z^  +  C 

=  r^^  +  2  [  +  %  ^  )  +  Tj  (y^  +  ^  n)  +  C  (Zg  +  i  C  )  ] 

2.2.3  ■ 


and  q  is  taken  to  be 

q  =  ^  +  n  (yg+^  n)  +  ^  2.2.4 

r 

e 

so  that  equation  2.2.3  becomes 
2 

(^)  -  1  +  2  q  2.2.5 

e 


The  first  term  on  the  right-hand  side  of  Equarlon  2,2.2  may 
be  expressed  as 


2.2.6 


where,  applying  Equation  •2.2.5, 

(^)  =  (1  +  2q)'^/2 


2.2.7 
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T^«n,  using  •  binomial  expansion,  Encke's  series  is  defined  by 


r  3 

i  - 


1  -  (1  +  2q) 

fq 


-3/2 


-■>  -  hi 

^  2.' 


^2,  3-5-7 


3.' 


2.2.8 


Finally,  the  substitution  of  Equations  2.2.6  and  2.2.8  into  Equation 
2.2.2  yields  Encke's  formula: 


^  (fqx-i)  +  C  2.2.9 

di; 

The  perturbative  term  x  typically  includes  the  second, 
third,  and  fourth  harmonics  of  the  earth's  gravitational  potential  and 
the  perturbations  due  to  the  sun  and  moon. 

The  procedure  employed  to  calculate  the  position  and  veloc¬ 
ity  on  the  reference  orbit  from  the  unified  parameters  follows  the  pat¬ 
tern  developed  in  Section  2.1. 
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2.3  DIFFERENTIAL  EXPRESSIONS 

This  section  deals  with  the  development  of  linear  differential 
expressions  relating  position  and  velocity  to  the  unified  orbit  para¬ 
meters  introduced  in  Section  2.1.  In  Section  2.4,  these  relationships 
will  be  extended  to  include  observation  residuals  referred  to  topocentric 
coordinates.  In  each  of  these  operations,  the  algebraic  manipulations 
are  tedious,  but  straightforward.  The  object  of  this  present  section 
is  vector  relationships  relating  position  and  veloci^  increments  Aj;  and 
Ar  to  increments  in  the  parameters  Ai,  and  A^-  In  order  to 
remove  redundancy,  the  latter  two  parameters  will  be  developed  in  terms 
of  related  quantities  which  recognize  the  unit  property  of  the  vector  !P 
and  the  orthogonality  of  ^  and 

Differential  expressions  provide  the  link  between  observation 
residuals  and  corrections  to  the  parameters.  The  development  starts 
with  the  vectors  of  position  and  velocity; 


r  =  X  P  +  -E  2 
■  ^P 


r  =  X  P  +-Ji’ 
■  -fp 


2.3.1 


2.3.2 


where,  for  convenience,  the  coefficients  of  will  be  denoted  as  and 
respectively.  Differentiation  leads  to: 

Ar  =  P  +  +  2.3.3 

Ai  =  P  Ax^ +2-^^  +  S  Aa  2.3.4 
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2.3.1  Differential  Expressions  for, Position  and  Velocity  Components  In 
the  Orbit  Plane 

In  order  to  facilitate  understanding  of  this  differential 
procedure,  those  derivatives  in  equations  2.3.3  and  2.3.4  which  pertain 
to  the  position  and, velocity  components  in  the  orbit  plane,  i.e., 

A^,  Ax^,  As,  will  be  derived  first  in  terms  of  corrections 
to  the  parameters.  Subsequently,  the  differentials  associated  with 
the  expressions  for  the  vectors  defining  the  orientation  of  the  orbit 
plane  will  be  discussed.  The  and  £  vectors  both  define  the  orienta¬ 
tion  and  serve  as  the  chosen  axes  of  the  orbit  plane. 

The  position  and  velocity  components  in  the  orbit  plane, 

^U)  »  yu)  ^  »  you  »  respectively,  are  defined. 


x^  =  q  - 

=  s:  -  i  U 
a 


2.3.5 


2.3.6 


•V  ^  - 

The  differential  expressions  for  x^  ,  S,  x^  and  S  are  ob¬ 
tained  from  the  differentiation  of  equations  2.3.5  and  2.3.6: 


^^00  =  ^ 

=  Ax  -  -  -  U  A  (i) 

=  f  [|Ar-A'S! 


2.3.7 

2.3.8 

2.3.9 

2.3.10 


where  the  various  terms  are  defined  below.  The  determination  of  AC 
and  follows  from  the  series: 


» 


14 


V 


ly  = 

'  31 

-xi-  + 

a-51  ^ 

1  *•  •  •  • 

•  / 1 

•Vi 

‘ic2 

'''it 
.  ^ 

C  = 

■  21 

a  *41 

I  1  *•  •  •  • 

•61 

Differentiating, 

we  obtain 

Au 

=  c 

AX'  -  ly 

a 

A(i) 

AX'  -  ?r 

a 

where 


3^5 

_  2 

3^7 

3 

'X9 

"  51 

a 

71 

+  "2 
a 

91  “  ‘ ' 

3j4 

^  2 

3(6 

3 

3(8 

*  41 

a 

61 

81  " 

From  the  formulae  for  e  and  r 


r  =  q  +  e  6', 
differentiation  yields 


Ae  =  -q  A  ( j)  -  J  A  q 

and 

Ar  =  Aq+eA'((  +  '<y  Ae 

or 

Ar  =  e^  AX'  -  (qfi'  +  e?!'^) 

A  (i)  +  (1 

15 


Finally,  the  ^  term  is  derived  from  the  expressions  for 


+  n(t-t) 
o  o 


fif.  q  S'  +  e  ly 

and  their  derivatives  which  follow 

A  jy  -  A  fif 


since  A  n  ■  0 


AM-qAX  +  eAU+XAq  +  UA 


or 


ASf--Afif  +  -(qty+eSf)  A(-)  -  -  A 


r  o  r 


a'  'a'  r 


At  thiSj^point  it  is  convenient  to  express  Eq.  2.3.7-10  in 
terms  of  A  fif  ,  A(— )  and  A  q; 

O  d 


Afif  A(-)  +  X  A  q 

^  m  o  aq  'a  q  ^ 


A'S'-^f  AM+?  A(-)  +  S  A 
m  o  aq  'a  q 


2.3.11 


2.3.12 


where 


X  •  -  S/r;  X  *-X  (eU  +  qU)+C;X  »1-X  S 
m  aqm'a^'  a  q  m 

.  (1  -  i  ?r)/r;  ^  ^  (e  +  q  Ij)  -  (S'  -  i  iT  ) 

m'  a‘^'’aq  m'  a  ^"^  '  a  a' 


ry  ‘V 

S  -  -S  S  2.3.13 
q  m 


and 


Ax„,»3[„  A\  A(i)  +'S  Aq 

*"  O  fl  Q 


O)  m 


■  'S  Afy  +  S'  A(— )  +  'I  Aq 

m  o  aq  'a'  ^  q  ^ 


2.3.14 


2.3.15 
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where 


'u  'J' ji  e  £ 

in  ”  3  r 

r 


q  3  r 

^  r 


2.3.16 


3  - - 2l  +eS)  S  =  -bb  -  r— 

in  '  q 


■^a,  =  i  ^  «»>  *  ^a)  -  «> 


To  more  directly  relate  the  A*S,  etc.,  to  the  chosen 

parameters,  the  quantity  Aq  Is  replaced  by  2  ^p.  Note  that: 

^•'^=P 


^ • aS^^ap 


Then,  from  the  expressions  relating  p,  q,  e  and  — 


q  =  p/{l  +  e) 


e  =  1  -  p/a 


one  obtains 


Aq  =  i  (JA  P)  +  A(i) 
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2.3.2  Differential  Relationships  Extended  to  the  Orientation  of  the 
Orbit  Plane 


The  orientation  of  the  orbit  plane  is  defined  by  P,  a  unit 
vector  directed  toward  perigee,  and  a  vector  normal  to  P  in  a  sense 
determined  by  the  direction  of  motion.  To  avoid  the  indeterminate 
behavior  of  Q  for  rectilinear  motion,  the  vector 

3:  =  -/p  2. 

has  been  considered  in  this  development. 

The  reduction  of  the  six  components  of  A  £  and  to  fewer 
quantities  leads  to  a  definition  of  £  in  terms  of  its  right  ascension 
a*  and  declination  6 *.  From  these  the  unit  vectors  A*  and  D*,  ortho¬ 
gonal  to  £  (see  Fig.  2)  are  calculated. 


P 


'  P  =  cos  6*  cos  a* 

.  P^  =  cos  5*  sin  ct* 

P  =  sin  6* 

\  ^ 


/ 


4*  < 


A^*  =  -  sin  a* 
A  *  =  cos  a* 

y 

A  *  =  0 
z 


t 


D*  J 


D^*  =  -  sin  5*  cos  cc* 

D^*  =  -  sin  5*  sin  0,* 

D  *  =  cos  6* 

2  ^ 
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The  differentiation  of  P,  A;*  and  D*  yields 

AP  =  A*  cos  6*  Aa*  +  ^  A6  *  2.3.23 

AA*  =  -P  cos  6*  Aa*  +  2*  sin  6*  Aa* 

AD*  =  -A*  sin  6  *  Aa*  -  I  A6  * 

The  angle  0*,  shown  In  Fig.  2,  allows  us  to  derive  the 
following  relationships. 

=  A*  cos  0*  +  D*  sin  0* 

W  =-A*  sin  0*  +  D*  cos  0* 

^  =  A*  c*  +  D*  s*  2.3.24 

where 

c*  =  ,/p  cos  0* 
s*  =■  Jp  sin  0* 

The  differential  expression  for  ^  becomes 

Aff  =  A*  A  c*  +  D*  A  s*  +  W*  Aa  *  +  p*  A6*  2.3,?5 

where 

W*  =  c*  (D*  sin  6*  -  P  cos  5*)  -  s*  A*  sin  5* 

P*  =  -P  s* 
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A  further  reduction  of  unknowns  in  the  expressions  for 
A  x^  ,  A^,  Ax^  and  A?  is  carried  out  by  noting  that 

1  •  2  -  i  Ap  -  c*  Ac*  +  »*A  »* 


2.3.26 


Incorporating  equation  2.3.26  into  equations  2.3.17-22  and 
substituting  these  along  with  equations  2.3.23  and  25  into  equations 
2.3.3  and  4  gives  the  followingj^ formulae  for^  A£  gnd  A£  ig  terms  of 
six  unknowns,  namely,  A?1  ,  A(“)»  Aa  »  A6  ,  Ac  and  As  . 


o  'a' 


[piC+tf&'l  ]  A(-)  2.3.27 

—  m  m  '  o  ap  ^  ap  'a' 


[i<  f\j  ^  f\j  "k  ^  .  "k 

c  (PX  +QS)+A  SjAc 
“  P  P  “ 

+  [s  (P  X  +  Q  S  )  +  D  S  ]  As 
“  P  P  “ 

+  [  A  x^  COS  6  +W  s]  ACX  +[D  +P  a]  A6 


Af=  [P^+'^^]Aa'  +  [pi  +^S  ]A(i)  2.3.28 

—  m  •^m  o  ap  *ap  'a' 


+  Lc  (PX+6s)+AS]Ac 
“  P  P  “ 

S  (PX  +QS)+D  sJAs 
-  P  ^  P 


a*+w*i  lAa*+  [d*  + 


k  k 

The  explicit  derivation  of  Ot  and  6  is  avoided  by 
expressing  the  involved  quantities  A*  and  D*  in  terms  of  the  given 
components  of  P: 


k 


* 

\  ■  0 


* 

D 


z 


-  (p'+  p 


2 

y 


* 

A 

X 


a  -P  /D 

y 


* 

z 


* 

D 

X 


P 

z 


* 

D 


y 


* 

A 


X 


P 

z 


1 


It  Is  obvious  from  Fig.  2  and  Eq.  2.3.24  that  c  and  s  may 
be  obtained  from  the  following  dot  products: 


^  •  A 


* 


* 

»  c 


-k  * 

2  •  D  as 

This  completes  the  derivation  of  the  differential  expressions.  In  the 
following  section,  the  topocentrlc  coordinates  of  the  observer  are 
Introduced,  and  scalar  differential  expressions  between  parameters  and 
residuals  In  the  observed  quantities  are  derived.  Details  of  the  structure 
of  the  IBM-709  program  are  reserved  to  Sect.  3. 


2.4  TOPOCENTRIC  RESIDUALS 

The  position  of  the  vehicle  referred  to  the  observer  Is 
related  to  the  positions  of  the  observer,  R,  and  of  the  vehicle,  £, 
both  referred  to  the  dynamical  center  by  the  fundamental  expression. 


£  •  i  +  R  ■  pL 


2.4.1 
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where  L  is  a  unit  vector.  (See  Fig.  3) 

The  differential  relationships  relating  residuals  in  the 
observed  position  to  improvements  of  the  orbital  elements  are  embodied 
in  the  derivative  of  equation  2.4.1, 


^  =  A  (pL)  «  ^  L  +  -  Ar  . 


2.4.2 


It  is  assumed  that  there  is  no  requirement  for  correcting  the  station 
vector,  R.  The  previous  section  developed  an  expression  for  in 
terms  of  corrections  to  the  orbital  parameters.  In  this  section  Ag 
is  translated  into  differential  expressions  involving  residuals  in  the 
observations.  In  the  discussion  to  follow,  let  the  subscript  c  signify 
that  the  particular  quantity  is  computed,  that  is,  obtained  from  the 
representation,  whereas  no  subscript  denotes  an  observed  quantity.  The 
residuals  will  be  taken  in  the  sense  "observed-minus-computed." 


The  observations  considered  in  this  study  are  p,  A  and  h, 
and  5,  and  p  .  The  differential  correction  program  processes  any 
combination  and  number  of  these  observed  quantities  taken  one  at  a 
time.  This  implies  that  although  a  complete  set  of  observations  may 
be  available,  for  example.  A,  h  and  p,  such  that  the  components  of  A^ 
could  be  obtained  directly,  the  more  general  method  employed  in  the 
differential  correction  program  considers  each  piece  of  information 
without  regard  to  the  others.  This  procedure  allows  Incomplete  data 
to  be  processed  and,  furthermore,  permits  Independent  weighting  or 
rejection  of  the  data  components. 


First,  if  only  one  or  two  of  the  a,  5,  and/or  p  are  available 
it  will  be  useful  to  define  the  auxiliary  unit  vectors  A  and  D  (Fig. 4 ) 
which  form  an  orthogonal  set  with  L.  D  is  the  north-pointing  tangent 
to  the  celestial  meridian  at  the  point  toward  which  L  is  directed.  A 
is  also  tangent  to  the  celestial  sphere  at  the  same  point  and  parallel 
to  the  equator  plane  so  that  it  completes  the  right-handed  orthogonal 
set  L,  A,  and  D.  The  components  of  L,  A,  and  D  are 

=  cos  6  cos  a  1 


L 


y 


cos  6  sin  a 


>  L 


2.4.3 


L  ■  sin  6, 
z  ' 
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North  Celestial  Pole 


Figure  3 

POSITION  RELATIONSHIP  OF  OBSERVER,  SATELLITE  AND  DYNMllCAL  CENTER 
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Celestial  Spher 
Centered  At 
Observer  / 


Figure  4 

VECTOR  RELATIONSHIPS 
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X  -sin  a 


Ay  =  cos  a  f  — 


2.4.4 


\  •  0, 


s  -sin  8  cos  a 


Dy  s  -sin  8  sin  Ot  L  D 


2.4.5 


D  m  COS  8, 


where  a  and  8  are  the  topocentric  right  ascension  and  declination  of  the 
object.  In  addition,  la  determined  from 


where 


^  “  £c  /  Pc 


£  »  r  +  R 

c  —  “ 


Pq  "  ^c  ^c 


2.4.6 


For  angular  observations,  the  residuals  in  the  direction 
cosines,  are  obtained  from 

41;  «  L  -  L 


Also  it  can  be  seen  from  Fig.  4  that 


2.4.7 


■  A  A  Otcos  8  +  D  A  8 


Successive  dot  products  of  A  L  with  A  and  D  will  lead  to  scalar  differ¬ 
ential  expressions  for  A  acos  8,  and  A&,  respectively. 


^  • 

A  a  A  OCcos  8 , 

2.4.8 

D  a  £,5 

2.4.9 

Similarly,  the  dot  product  of  ^  with  L  relates  a  slant-range  residual 
to  improvements  of  the  orbit  parameters. 

•  L  ■  2.4.10 


Differential  expressions  involving  observations  in  a  horizon- 
oriented  coordinate  system  may  be  obtained  in  a  similar  manner  by 
defining  the  horizon  system  unit  vector  triad  L,  A,  and  D,  whose 
components  are: 


«-cos  A  cos  h 

s  sin  A  cos  h  ^ 

L  ,  *  sin  h, 

T.h  '  y 


L 


A  ,  a  sin  A 
xh 


'V 

Ayh  =  cos  A 


2.4.11 


2.4.12 


and 
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=  cos  A  sin  h 


=  -sin  A  sin  h 


cos  h, 


2.4.13 


where  h  is  the  altitude  or  elevation  angle  and  A  is  the  azimuth, 
measured  in  the  positive  sense  to  the  east  from  north.  These  horizon 
system  components  are  rotated  into  the  equatorial  system,  to  which 
the  components  of  are  referred  by 


L  ■L.SjL.Ej.L.Z  ''i 
X  xh  X  ^  yh  X  ^  zh  x 


=  A,S+A.E-'  A.Z 
X  xh  X  yh  X  zh  X 


(x — ►y,  z.)  2.4.14 


D  sD.S+D.E+Du^^ 
X  xh  X  yh  X  zh  x' 


The  components  of  the  auxiliary  unit  vectors  E  and  Z,  directed  to 
the  south,  east,  and  zenith,  respectively,  are 


r  sin  0  cos  0 


S  ■  sin  0  sin  0  I  S 

y  i  - 


2.4.15 


S  ■  -cos  0  , 

z  ' 


E  s  -sin  0 

X 


E  =  cos  0 

y 


E_  -  0, 


2.4.16 
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and 


Z 

=  cos  <t>  cos  6  ] 

1 

X 

Z 

y 

] 

=  cos  sin  Q 

.  z 

2.4.17 

z 

■  sin  'P  , 

z 

where 

9  is  the  observer's  astronomical 

latitude  and 

6  is  the  local  hour 

angle 

of  the  vernal 

equinox.  Then 

■  "X  A  A  cos  h  +  Ah. 

2.4.18 

Scalar  differential  expressions  for  A  A  cf^s  h,  ^nd  Ah  can  now 
be  obtained  by  forming  the  dot  products  of  with  A  and  D,  respectively: 


‘  A  •  AA  cos  h 


2.4.19 


» 

and 


•  D  •  Ah. 


2.4.20 


Thus,  residuals  in  the  observations  of  angular  position  in  the  horizon 
coordinate  system  can  be  related  to  corrections  to  the  orbit  parameters 
through 

Range-rate  residuals  can  be  related  to  the  parameters 

through 


£  •  Aq  =  P  A  p 

2.4.21 

£  •Z^+£  •  -  pAp  +  p 

A  p 

2.4.22 

or 

pAp-£-Z^  +  [£  -  £]• 

A£ 

2.4.23 

The  corrections  to  the  orbit  parameters  are  introduced  through  ^  and 
as  developed  in  Sect.  2.3. 
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SECTION 


3 


EXPERIMENTATION 


The  unified  ephemeris  and  differential  correction  theory 
has  been  subjected  to  experimental  verification  to  determine  the 
V '.curacy  and  rate  of  convergence  with  various  types  of  simulate!  observations. 

EPHEMERIS  PROGRAM 


The  numerical  integration  program  was  checked  against  the 
Encke  Moonshot  program.  The  major  differences  between  the  two  methods 
are  the  choice  of  parameters  used  to  define  the  two-body  reference 
orbit  and  the  argument  used  for  the  numerical  integration,  i.e.,  time 
for  the  unified  Encke  program  and  eccentric  anomaly,  E,  for  the  Encke 
Moonshot  program.  The  results  of  the  comparison  are  listed  in  Table 
1. 


Further  experimentation  is  being  conducted  to  determine 
integration  step-size,  time  for  integration,  etc. 

DIFFERENTIAL  CORRECTION  PROGRAM 

The  program  has  been  verified  by  correcting  a  hyperbolic 
orbit  to  one  that  is  elliptical  with  various  combinations  of  range 
and  angle  data.  Tables  2  through  4  give  the  sizes  of  the  r.m.s. 
residuals  and  list  the  values  of  the  parameters  M  ,  i,  P  and  ^  after 
each  of  four  corrections.  The  results  of  Table  2  were  based  on  32 
range  observations  spaced  over  7  hours  and  18  minutes.  Over  the 
same  time-span.  Table  3  involves  15  observations  of  ^'and  6  and  17 
observations  of  A  and  h,  and  the  same  number  of  angular  observations 
were  used  with  32  range  observations  in  Table  4.  Table  5  is  a 
reproduction  of  the  program  listing  of  the  residuals  in  the  range  and 
angular  observations  after  one  correction. 
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Further  verification  of  the  program  was  carried  out  by 
correcting  from  an  elliptic  to  a  parabolic  orbit  and  from  a  hyperbolic 
to  a  parabolic  orbit.  The  success  of  this  experimentation,  orbit 
correction  through  unit  eccentricity,  satisfies  the  definition  of  the 
unified  theory.  Tables  6,  7  and  8  list  the  values  of  the  parameters 
for  the  elliptic  case  and  the  r.m.s.  residuals  in  the  observations 
after  each  correction  for  both  the  elliptic  and  hyperbolic  runs. 

In  all  cases  the  observations  were  obtained  over  a  6  hour  and  42  min. 
period.  Thirty- two  range  observations  were  used  in  Table  6,  13  sets 
of  Otand  6  and  15  sets  of  A  and  h  for  Table  7,  and  a  combination  of 
32  range,  13  a  and  5  and  14  A  and  h  observations  are  included  in  Table 
8.  Tables  9  through  12  demonstrate  the  use  of  the  program  with  range 
rate  simulated  observations. 


COMPARISON  OF  ENCKE  MOONSHOT  PROGRAM 
TO  UNIFIED  ENCKE  INTEGRATION  PROGRAM 
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At  implies  the  interval  of  time  since  start  of  integration. 


RANGE  OBSERVATIONS 
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